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In this work, the Friedman equations for hadronic matter in the Robertson-Walker metric in the 
early Universe are obtained. We consider the hadronic phase, formed after the hadronization of the 
quark-gluon plasma, that means times from 10 _6 s to Is. The set of equations will be derived and 
the behavior of the system will be studied considering one approximate analytical solution. 


In the past few years, at CERN [H,@ and at RHIC 
i-1 very important results have been obtained. It has 
been observed that when matter is submitted to the ex¬ 
treme conditions of the high-energy heavy-ion collisions 
(high temperatures, T > 150 MeV, high energy densities, 
exceding 1 GeV/frn, and small barionic chemical poten¬ 
tial p, ~ 0), it is believed that the quark-gluon plasma 
phase is produced. At the theoretical level this under¬ 
standing is improving @|-[1], and gradually, the behavior 
of the QCD phase diagram is being clarified. 

These conditions are believed to be close to the ones 
found in the early Universe, so, we may be tempted to 
try to understand the early Universe in a similar way [§!]- 
11]. In high-energy collisions, after the QGP creation, 
it expands and cools, and the hadronization occurs dur¬ 
ing this expansion, determining the freeze-out process 

mM- 

The hadronic phase produced in high energy heavy- 
ion collisions has been successfully studied in [1J]-[16|, in 
terms of effective relativistic nuclear models considering a 
mean-field theory by means of a Walecka model 0 , and 
its main features could be explained with this procedure. 
This fact support the idea of using this model, at least 
as a first approximation, in this work. 

On the other hand, one must consider hot hadronic 
matter in a curved space-time. In 0,0 this kind of 
system has been considered, but for studying the neutron 
stars structure (T ~0). 

So, in this work, we investigate the evolution of the 
early Universe, by the Friedman equations, considering 
an effective relativistic mean-field model for the hadronic 
phase, with the strong interactions given by the exchange 
of <7, w and p mesons. We suppose that this model is valid 
after the hadronization and before the nuclei formation, 
that means times from 10 _6 s up to Is, approximately. 
The equations will be obtained in this framework and 
one approximate analytical solution will be shown. 

To study the evolution of the Universe we will consider 
the Robertson-Walker metric 


ds 2 = dt 2 — R 2 ( t) 


dr 2 


1 — kr 2 


+ r 2 dd 2 + r 2 sin 2 9 dtp 2 


(1) 

that assumes an isotropic and homogeneous hadronic 


phase. 

In this phase, we will consider the a — oj — p model 
adapted to the curved space-time in a way similar to the 
one shown in fl9| . but now, considering the Robertson- 
Walker metric ©. The Lagrangian density is 

£ = £& + £ e + C v + £ CT + Cuj + Cp + Cj + C-int + £g , (2) 

where the first term, Cb, describes the baryons (protons 
and neutrons) and, in the same way, the terms C e and £„ 
are related to the electrons and neutrinos respectively. 
The term refers to the Lagrangian densitiy of the 
free-scalar boson field <r, which describes the long range 
attractive interaction. The terms C u and C p are the La- 
grangians of the free massive vector boson fields w and 
p, that describe the short range interactions. Now, the 
term £ 7 is the Lagrangian of the electromagnetic field. 
Finally, Cmt represents the interacting part of the La¬ 
grangian density and the term C g is the Lagrangian of 
the free gravitational field. The quantities g a , < 72 , < 73 , <7^ 
and g p are the coupling constants, and can be found in 
the in the references [23]-[111, 03 is the Pauli matrix, and 
e is the proton charge. Other observation is that is 
the covariant derivative of the respective object that it 
acts. These Lagrangian densities may be expressed as 


Cb = * (Vv M - Af 6 ) T, (3) 

C e = Ip e - m e ) ip e , (4) 

£<x = iv^crV^CT - im 2 cr 2 — U (a ), (5) 

C u = ( 6 ) 

C p = -\~ p „~^ + l -m 2 pPp p^ ( 7 ) 

= ( 8 ) 


C-int = -gu, i &'f‘u l Pb - -pp’Fq'Va/cv’P 

^ + eV> e 7 M AAe, (9) 
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with 


M b 




1 — [IIP 



T 

U(a) 


( rn p - g a a 

\ o 


0 ) 

m n - g a cr J 


( 11 ) 


(rn p 0 \ 

V 0 fh n ) ' 


d p u}„ - d u u} p , 

dvPfl: 

^pA^ ^l^A p , 



p —> proton, n —> neutron 



+ 4^ a 
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( 12 ) 

( 13 ) 

( 14 ) 

( 15 ) 

( 16 ) 
( 17 ) 


Two observations must be made about this formula¬ 
tion: we will consider the difference of mass between the 
proton and the neutron and we will not show the neu¬ 
trino’s Lagrangian. The reason to consider the difference 
of mass is because it causes a discrepancy between the 
neutron and proton densities that we desire to capture 
with the model. Now, the reason for not showing the 
Lagrangian of the neutrinos is just because we will not 
be concerned with a model for neutrinos in a first ap¬ 
proach, this study deals with hadronic interactions in the 
hadronic phase and we can approximate its contribution 
as an ultra-relativistic Fermi gas as it will be made be¬ 
low. The spin connections that appears in the fermion’s 
covariant derivatives, that should be considered in this 
kind of picture, will be neglected in this formulation as 
its effect is negligible in a mean field theory, and this 
study will be left as a refinement for future works. 

Now if we define 


($*)=p a , (18) 

<*7^> = K = + JZ, (19) 

(WyVsvk) = Jg = (20) 

(*7 M *) = J p > ( 21 ) 

(Te^e) = j (22) 

(a) = <j>, (23) 

(n # “') = V , (w^) = V», (24) 

(S » V )=B' IV , (25) 

{Fn = , (Ap) = A p , (26) 


we may derive the equations of motion by applying the 
mean-field approximation 


(V=9d»<t>) = g<f,p s ~ m\<$> - ( 27 ) 


1 


1 


dp {V^gV^) = g v r b - m 2 v V v 

v 2 t rv 

= gyribU — m v V , 

i 


(V=gB^) = - gB (j; - r n ) - m 2 B B - 


= -jPB (rip - n n ) u v - m 2 B B v ,( 29 ) 

-±=dp W=gF^) = e (j; - r e ) 

= e(n p — n e )u v , ( 30 ) 


where is the average quadri-velocity of the charge car¬ 
riers, and the densities are 


ps = J2 


(2nY 


d 3 k 


mi 


Vk 2 


(fi+ + fi-) ) ( 31 ) 


n b 


= Y 77T73 f d3k (^+ “ /i-) = n p + n n, (32) 
-ki J J 

f d 3 k(f e+ -f e _), ( 33 ) 


p,n 
2 


(2tt) j 


where 


f izt e /3 (EiTvi) + i ’ 

f 1 
Je± e P(E eT v e ) + X ’ 


i = p,n 


( 34 ) 

( 35 ) 


that are the Fermi-Dirac distributions with index “+” 
for particles and ” for anti-particles, v,-, and v e are the 
effective chemical potentials that have the form 


9vV ° ~ 2 9bB °’ 

( 36 ) 

gvV° + \g B B°, 

( 37 ) 


( 38 ) 


in a homogeneous and locally neutral system (for electric 
charge) and with u M = (1,0) (that is the case that we 
have here by virtue of the constraint imposed by the met¬ 
ric). Note that the factor “2” in each density is the spin 
degeneracy factor and we also have changed the index in 
the equations (|271) - (I5U1) for correspondence with the field 
which it refers. The equations derived above show how we 
can study a medium, composed by interacting hadrons, 
considering the mean-field approximation. In this frame¬ 
work, the metric is also considered, and the composition 
of this phase also has effect on the determination of the 
metric. 

Now, we can apply these ideas in order to describe the 
metric in the hadronic phase of the Universe, that means 
times above 10~ 6 s and below Is. For this purpose, as a 
first approximation, we will consider a isotropic and ho¬ 
mogeneous Universe, described by the Robertson-Walker 
metric ©, that determines the local neutrality for elec¬ 
tric charge. The constraints lead us from the equations 
m-m to the set of equations 


d 2 (f) 

dt 2 


, R dcj) 
R dt 


dU 

d(j> 


— g<t>Psi 


myV° = g v J b = g v n b , 
m 2 B B° = g B J$ 


— gB (rip Tin) i 


( 39 ) 

( 40 ) 

( 41 ) 


( 28 ) 
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where we suppressed the electromagnetic field equation 
because it vanishes if homogeneity and isotropy are im¬ 
posed. Another important condition is the conservation 
of the baryon number that is expressed by = 0 

and gives the result 

n b = n b0 . ( 42 ) 


and the definitions e b +eB+e r = £t and p b +PB + ^f = Pt 
to obtain 


/ der\ d/3 / <9et\ dz 

V <9/3 y ^ dt + V dz J ^ dt 



+ 9<f>Pa 


dcj) 

dt 


—3— ( pt + £t) ,(46) 


where n b o is the initial baryon density (when the hadron 
epoch begins). 

Supposing that the matter may be considered as a per¬ 
fect fluid, we use an energy-momentum tensor of the form 
T^ v = (e + p) u^u v — pg^ with u^ = (1, 0 ). Making use 
of the Lagrangian densities determined above, we may 
obtain 


e — £,/> + £y + e.B + £e + £j/ + £7 + £b 

= 1 {T) +r ^ 2 + £/w 


1 9v 


2 rriy 


Rc 


1 9 b 


43 


2rriA 


ZV Z I 'd I I c7 1 v Zrr r 

0^.2 n b0 ( ^> ) 4" o ™2 ( n p n ‘ n ) 4" ion 1 ^ 


2 


120 


+ 77^3 H [ d 3 kyjk 2 + m 2 (/i++/i_) : 

( 2 ' 7r ) d 

P = P</> + Pv + PB + Pe + Pv + P-t + p b 

= l(¥r) - l 171 ^ 2 - u W) 


(43) 


dt 


1 ffy 


2 TOy 

1 2 


i?r 


1 


43 


2/ti4 


^ V Z / u 1 X vts ( 1 2 /t>< 

o n b0 ( ^> J + o™2 V n P n n) + q^h ^ ^ 


2 


3 ( 27 t ) 


3 

p,n * 


d J /c- 


\//c 2 + mf 


360 

(/*++/»-) • (44) 


But is necessary to make some observations. As we said 
before, the electromagnetic field vanish if we use rigor¬ 
ously the Robertson-Walker metric. However the contri¬ 
bution of the photons is significant in this phase of the 
Universe. Then we included the photons as a boson gas, 
that is a good approximation, as electromagnetic inter¬ 
actions plays a role of small corrections in the hadronic 
phase (fact that is widely used in the study of high-energy 
heavy-ion collisions). For similar reasons, the electrons 
and neutrinos are approximated by a Fermi gas without 
chemical potential. These approximation lead us to con¬ 
sider p p = p n = p b in the /3-equilibrium. To conclude, 
these are some approximations that we made, for simplic¬ 
ity, in order to sketch the model and obtain results, and 
may be easely included inside this formalism in a future 
work, just by adding the correspondent terms in the La¬ 
grangian density. These three components are expressed 
together in the term proportional to T 4 . 

With the energy density and pressure expressed as 
above, we can use the equation 

= -3pR 2 (45) 


d (eR 3 ) 


where z = and /3 = ^. Note that the index in 
derivatives are the variables that are left constant. 

Now we can complete the set of equations with the 
Friedman equation (the first equation written below) and 
get 


d 2 (f> 77 dcj) 

dt 2 77 dt 

der \ d/3 

<9/3 Jzj dt + 


n 2 



dU 
dcj) 
dz 
g,(/> dt 



9<t>Ps 


dcj) 

dt 


n b 



(47) 


9<j>Ps, (48) 


-3— ( pt + £t) ,(49) 


n b o 

W ’ 


(50) 


that are the equations for the hadronic phase of the Uni¬ 
verse, the main objective of this work. The numerical 
solution of these equations will be published soon. Here 
we have used the fact that the constant fc, which should 
appear in the Friedman equation, can be taken as 0, what 
leads to 77. = This result will dispense us from the 
necessity of discovering the value of i?o, which demands 
the knowledge about the previous phases of the Universe. 

Although the complexity of the equations, we can get 
analytical results if we use some approximations in order 
to study the main characteristics of the solutions. First 
we consider only the term with T 4 in the energy density 
and pressure (as we are in the radiation-dominated Uni¬ 
verse) and, through the equation (l45l) . we get T = ^ 
and 


77 = -\/l2,3 • 10 3 t + 1 (51) 


where the time is given in seconds. In a non-relativistic 
approximation 


E i ~m i + 


1 k 2 

2 rrij ’ 


(52) 


Pi = rrii+p^ R , and considering a very small contribution 
of the (j) field (and of the vector fields), we get 


fi+ 


fi- 


r ( 1 ±2-—,,*' --. 
e P{2 mi V* ) + 1 


,NR~\ 

1 


i = p,n. 


(53) 


dR 


/(l^+A WR +2 m t ) + 1 


(54) 
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Now, with the approximation 


r xX ~ ldx we, 

/ -1 x , 1 - r ( A )U 

J o Z 1 e x + 1 

(55) 

where T (A) is the gamma function, we have as 
particles and anti-particles 

density of 

2 (2-77771,) 5 3 

n l+ = 3 zT 2 e t , 

(2tt) 

(56) 

2(2tt rmY -lrp- -?x. 

(27r) 3 

(57) 

So we can write 


2 (27TTOj) 5 , _ n 3 

rii = -,— (z-z ) T 2 e f. 

(2n) 3 V ’ 

(58) 

With the condition rib = n n +n p = = n b o ( 

^) 3 and 


noting that z — z 1 = 2 sinh (/3/i fc ), we get 

Pb = Tarcsinh (a), (59) 


where 


a = 


{2nfn b0 (jr) 


e t 


4 (27r m n T ) 2 ( 1 + 


(s)' 




(60) 


Then we have the densities given by 

n l+ = 2 ( 2 J W ^) 2 T i e -^ +arcsinh(a) , (61) 


rii- = 


(Inf 

3 

2 (27rrrii) 2 ^ 1 a^_^j._ arcsinh ( a ) 

( 2 7r ) 3 


(62) 


The graphic of these expressions is shown in Fig. 1. As 
we can see, at high temperatures, the densities of the 
considered particles are very similar, but as the temper¬ 
ature decreases, p and n decreases, and only protons and 
neutrons remain. The densities of p and n, that at high 
temperatures are the same, present different values as 
the temperature becomes small. All these characterisitics 
appears in the results of 11], where the calculations are 
made with a different model. So, we may conclude, that 
even with the approximations made, the final results are 
consistent. 

In this work we have studied the first second of the 
Universe (after 10 -6 s), that is, the hadronic phase. We 
obtained the analog of the Friedman set of equations, 
considering hadronic matter described by a relativistic 
mean-field model in the Robertson-Walker metric. This 
model opens possibilities for many variations and refine¬ 
ments, that can be made in a straightforward way. In 
order to check the theory, a non-relativistic approxima¬ 
tion has been made and consistent results have been ob¬ 
tained. The main characteristics for the p, n, p and n 
have been described. 



FIG. 1. Baryons and anti-baryons densities. The dashed 
(with smaller dash) line is the proton density, the dashed 
interspersed with points line is the neutron density, the con¬ 
tinuous line is the anti-neutron density and the dashed (with 
biggest dash) line is the anti-proton density. 
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